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Simple analytic considerations are applied to recently discovered patterns in a generalized Fisher 
equation for population dynamics. The generalization consists of the inclusion of non-local compe- 
tition interactions among individuals. We first show how stability arguments yield a condition for 
pattern formation involving the ratio of the pattern wavelength and the effective diffusion length 
of the individuals. We develop a mode-mode coupling analysis which might be useful in shedding 
some light on the observed formation of small-amplitude versus large-amplitude patterns. 

PACS numbers: 87.17.Aa, 87.17.Ee, 87.18. Hf 
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We have shown recently that the Fisher equation used 
frequently for investigations of biological or ecological 
systems, when generalized to include spatially non-local 
competition interactions, leads to interesting patterns in 
the steady state density 0, 0] ■ In this Note we attempt 
to shed some analytic light on the formation of these pat- 
terns. The original Fisher equation 0,0] is 



du (x, t) 
dt 



DV 2 u (x, t) + au (x, t) - bu 2 (x, t) , (1) 



where u (x, t) is the population density of individuals un- 
der investigation (bacteria, rodents, etc.) at position x 
and time t, and D, a, b are, respectively, the diffusion 
coefficient, population growth rate, and competition pa- 
rameter. The generalized equation 0,0 features compe- 
tition interactions linking u{x,t) at point x with u(y,t) 
at point y through an influence function f a (x, y) of range 



(T, 



du (x, t) 
dt 



= DW 2 u(x,t) +au(x,t) 

-bu(x,t) / u(y,t)f a (x,y)dy, 
Jn 



(2) 



f2 being the domain for the non-local interaction. 

In 0], we found that the introduction of the finite- 
range competition interactions gives rise to the emer- 
gence of patterns in the steady state density u{x) with 
the following features: 

• No patterns appear 0] in the two extremes of zero 
range (in which the generalization reverts to the 
Fisher equation) and full range (in which the pop- 
ulation density is linked equally to all points in the 
domain) . 



• The pattern structure depends crucially on features 
of the influence function, specifically, its cut-off 
length and its width. 

• Even when patterns appear, their steady-state am- 
plitude can change abruptly from substantial to 
negligible as the parameters of the system are var- 
ied. 

• The critical quantity determining the separation 
of large-amplitude patterns from small-amplitude 
ones appears to be the ratio of the cut-off length of 
the influence function to its width. 

These findings raise two questions. Why do the pat- 
terns form at all? And, what causes the separation of the 
large-amplitude patterns from the small-amplitude ones? 
Both questions are interesting. The first is amenable to 
understanding via standard stability analysis consider- 
ations. The second is more difficult but might be ap- 
proachable through a mode-mode coupling analysis, as 
we show below. 

In order to address the first question, consider a 1- 
dimensional version of Eq. J3J f° r simplicity, and substi- 
tute in it 



u(x, t) = uq + e cos(kx) exp((pt). 



(3) 



Here uq is the homogeneous steady-state solution a/b. 
Considering periodic boundary conditions, and retaining 
only first order terms in e, we obtain the following dis- 
persion relation between the wavenumber k of any mode 
of the pattern and the rate ip at which it tends to grow: 



-Dk 2 - aT(k). 



(4) 



In this expression, the influence function (assumed to 
be even) is represented by its cosine (Fourier) transform 
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!F(k) defined as 



J-{k) = I cos(kz)f a (z)dz. 
Jn 

Stable steady-state patterns require that 



A > 2tt, 



D 



(5) 



(6) 



where A = 2ir/k is the wavelength associated with the 
k— mode of the Fourier expansion of the pattern. 

Condition allows us to check for the existence or 
absence of inhomogeneity, i.e., patterns, in the steady 
state. We see from © that the Fourier transform of the 
influence function at the wavelength under consideration 
should be negative for the patterns to appear and that 
its magnitude should be large enough. One way of un- 
derstanding this condition is to recast it as requiring that 
2ir times the 'effective diffusion length' should be smaller 
than the wavelength for the patterns to occur. By the ef- 
fective diffusion constant is meant D divided by —F(X), 
which is a factor decided by the influence function, and 
by the diffusion length is meant the distance traversed 
diffusively in a time interval of the order of the inverse of 
the growth rate. If the influence function is smooth such 
as in the case of a Gaussian in an infinite domain, the 
Fourier transform is positive and no patterns appear. A 
cut-off in the influence function produces oscillations in 
the Fourier transform which can go negative for certain 
wavelengths. The reported finding that the cut-off na- 
ture is essential to pattern formation can be understood 
naturally in this way. 

Let us consider, in turn, three cases of the influence 
function which we have used in our earlier investigations 
PJ: square, cut-off Gaussian, and intermediate. 

First we take 

f{x -y) = ^{9[w - (x - y)]9[w + (x - y)}}, (7) 

where 9 is the Heaviside function. The influence function 
is thus a square of cut-off range measured by w from its 
center. We will consider the case here when the range w is 
smaller than, or equal to, the domain length L. Equation 
<@J then involves an integral from to w, and gives 



9 



_ aS m(kw)_ Dk2 
kw 



(8) 



In terms of dimensionless parameters 
k' 



kyjD/c 



V 



= a/D, 



we have 



, = sin(fc / 77) _ fc/2 
k'rj 



(9) 
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FIG. 1: The dispersion relation 10 between the dimensionless 
growth exponent <p and wavenumber k plotted for different 
values of the ratio n of the influence function range to the 
diffusion length (see text). Values of r\ are 50 (solid line), 10 
(dashed line), and 2 (dotted line). Patterns appear for those 
values of k! for which <p is positive. 



which we plot in Fig. 1 for three different values (50, 
10 and 2) of the ratio rj of the width to the diffusion 
length (not effective diffusion length). For the third case 
there are no patterns: diffusion is strong enough to wash 
them out. For the intermediate case, patterns can occur 
with wavelengths corresponding to values around k' f» 
0.4 while for the 77 = 50 case, they occur around kl w 0.1. 

The earlier finding Q, Q that no patterns appear for 
extremes of the range of the influence function is clear 
from Eq. ©. As the influence width vanishes, i.e., as 77 
goes to zero, both terms in tp are negative and there can 
be no steady state patterns: we recover the solution for 
the local limit, corresponding to Eq. QJ, when w — > 0. 
Since the boundary conditions are periodic in a domain 
of length L, there are only the allowed values k = rnr/L 
of the wavenumber. Therefore, in the opposite limit of 
full range, i.e., w — > L, the sine term vanishes, <p' = — fc' 2 , 
and again there are no patterns. 

Precisely the same qualitative behavior occurs for 
other non-square influence functions such as the Gaus- 
sian with a cut-off, i.e., for 



f(x -y) = 



1 



/7rerf (w /a) 



exp 



y 



(10) 



We again consider the case when the cut-off length does 
not exceed the domain length. This leads to the Fourier 
transform of the influence function involving an integral 
from to w. In these as well as other cases considered, 
it should be appreciated that the domain length L, if 
taken to be smaller than the cut-off length, becomes itself 
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the cut-off length: factors such as kw appearing in the 
Fourier transform become then kL instead. 

For this cut-off Gaussian case, the square case disper- 
sion relation © is replaced by 

aexp[-(fco-/2) 2 ] 



erf 



2 erf (id/it) 
w 

~u 2" J 



ika\ 



erf 



ika 



Dk 2 



The dimensionless version is replaced by 

cxp(-fc' 2 /3 2 ) 



2erf(a) 



{erf(a - ik'0) + erf(a + ik'(3)}-k K 



(11) 

Here a and 2(3 are the ratios of the cut-off length to the 
range and of the range to the diffusion length respec- 
tively: 



w/a, 



2(3 = (TyQD. 

What is analogous to r] in the square case is their product 
2a(3 — w^/a/D. Plots which are essentially the same as 
those in Fig. 1 can be drawn for this Gaussian case. 

It is interesting to note that, while there is a single 
quantity 77 in the square case, there are two quantities, a 
and (3, in the cut-off Gaussian case. This arises from the 
fact that, although there are generally two lengths asso- 
ciated with any influence function, the cut-off length and 
the width, the latter is infinite for the square case. The 
width has been defined in Ref. Q as being inversely pro- 
portional to the second derivative of the influence func- 
tion evaluated at its central point, and has been denoted 
by the symbol X. The cut-off length measures the dis- 
tance beyond which the influence function is exactly zero 
and has been denoted 0] by £ G . For the cut-off Gaussian, 
this £ c = x c = w. The symbol x c has been used in Ref. 
PJ and w in the present Note. The width S obeys X = a 
for the Gaussian case and X = 00 for the square case. 

Both the cut-off Gaussian and the square can be ob- 
tained as particular cases of a general function 0, 
which we call the intermediate influence function: 



/0 - v) 



r(l/r + 3/2) r(x-yf 
y/nwT(l/r+l) (2 

{6[w - (x - y)]0[w + (x - 



- 3r)a 2 

y)}}- 



l/r 



(12) 



For r — ► the intermediate influence function goes to 
the cut-off Gaussian, while for r — * 00 it yields he square 
function. The cut-off length of the influence function is 
given by 



w 



(2 + 3r) 



(13) 




FIG. 2: The dispersion relation ll'it between the dimension- 
less growth exponent ip and wavenumber k plotted for the 
intermediate influence function. Values of rj are as in Fig. 1: 
50 (solid line), 10 (dashed line), and 2 (dotted line). Patterns 
appear for those values of k' for which tp is positive. 



We will follow the notation 

v = l/r + 1/2, 

and evaluate the Fourier transform of the influence func- 
tion by calculating the integral 0] 



2w- 2y T(v+ 1) 



cos 



(ks) [1 



ds 



— r(u + i)j v (kw), 

kw , 



(14) 



for 



k > 0, w > 0, Re[v] > -, 

r and J being the gamma and the Bessel functions re- 
spectively. The dimensionless dispersion relation analo- 
gous to (2J is, for this general case, 
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r(^ + l)J,(fc'r,)-fc' ; 

c 1] J 



(15) 



Here, as in the square case, 77 = w^/a/D is the ratio of 
the influence function width to the diffusion length. 

It is straightforward to obtain the two limits, square 
and Gaussian, from this dispersion result 115t for the 
intermediate influence function. In Fig. 2 we plot the 
intermediate case for v = 1 and see the same general 
behavior as in the Gaussian and the square counterparts 
(see, e.g., Fig. 1). Steady state patterns appear only 
around kl = 0.1. 
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Having understood the appearance of the patterns, we 
now come to the second issue mentioned in the intro- 
duction: the transition from small-amplitude to large- 
amplitude patterns 0]. This is much more difficult to 
address analytically. We present here a procedure that 
we believe has the potential to shed some light on this is- 
sue. We substitute the Fourier mode expansion of u(x, t), 



u (x, t) 



^ A n (t) cos(k n x), 



(16) 



in explicitly noting that k n = irn/L, and using the 
orthogonality properties of trigonometric functions, ob- 
tain separate equations for the n = mode, 



dA 
dt 



aA 



bAl 



n=l 



A; 



■J" (kn 



and for other modes n ^ 0: 



dAn 
dt 



- Dk n A n -\- (iA r 



(17) 



(18) 



-bA An [l + T{k n )\ 
j -T(k 3 ) 



A A 



3 = 1 

°° A. A. 

^ 2 

j=n+l 

[Hk^+Hkj-n)]- 

Equations ijTTjl and itlSll are the complete set of equations 
for the evolution of the amplitudes of all modes in the 
non-local problem given by Eq. (01 . The appearance of 
stable patterns only for those values of k n for which ip is 
positive as seen in our Figs. 1 and 2, suggests that we 
envisage an interaction between only two modes, the zero 
mode and the one whose growth we examine, say n = m. 
In a situation as in the plots shown in which ip > only 
for a small k— range, the discrete nature of the allowed k 
values could lead to only a single non-zero mode lying in 
the stable range. Then we would have only two coupled 
nonlinear equations for the mode amplitudes, 

J A A2 

^ = a A -bAl-b^(k m ), 



dA r 



dt 



= -Dk 2 n A m + aA m - bA A m [1 + T(k m )} 



which lead, in the steady state, to 

a - Dk?„ 



A 



b [1 + F[k v 



A^y. 2 



bAl - aA 
bT(k m ) 



(19) 
(20) 



Substitution of the zero mode amplitude in A ri 
an explicit expression for the latter: 



gives 



A 



T(k m )b 



Dk 2 m ) 2 a(a~Dk 2 J 



b{l+T(k m )f b{l+T{k m ) 



(21) 

Preliminary investigations lead us to believe that the 
appearance of small and large amplitude patterns for 
different parameter regimes might emerge from consid- 
erations of (1211 . The two-mode approach is, however, 
plagued by the fact that variation of the parameters of 
the system can make the approach invalid in one regime 
even if it is valid in another. Surely steady state pattern 
amplitudes can always be obtained by a simultaneous so- 
lution of the algebraic equations obtained by putting the 
left side of Ill7fl and II18H to zero. This is a straightforward 
numerical task if we can make the reasonable assumption 
that the number of modes to be considered has a finite 
cut-off. Such a cut-off is an obvious consequence of the 
fact the Fourier transform of a typical influence function 
disappears for high values of n. These and related studies 
will be reported in the future. 
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